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Kinetic analysis of the olive pomace thermal degradation in the temperature range of interest for torrefaction 
was performed by using non-isothermal thermogravimetric measurements at different heating rates, ranging 
from 2 to 40 °C/min. A comparison is presented between two selected integral isoconversional methods, 
i.e., the nonlinear Vyazovkin incremental approach, which is more accurate but time-consuming, and the linear 
Ozawa-Flynn-Wall (OFW) method, which is less accurate but computationally simpler. Results show that the 
values of the activation energy by the OFW method are consistent with the ones provided by the Vyazovkin ap¬ 
proach. This implies that the OFW method, more user-friendly compared to the Vyazovkin procedure, is suitable 
for studying the torrefaction kinetics of residual biomass, such as olive pomace. The reliability of the OFW method 
was further confirmed by the successful application of the derived kinetic data to reproduce (i.e., predict) exper¬ 
imental TG curves not included in the kinetic computations. 

© 2014 Elsevier B.V. All rights reserved. 


1. Introduction 

In the view of shortage of fossil fuels and with the increasing concerns 
regarding human impacts on the environmental, renewable energy 
sources and waste materials play an important role as a viable alternative 
to fossil fuels for both energy generation and production of chemicals [1 ]. 
Now on the brink of commercialization, torrefaction - a relatively 
new thermochemical treatment of biomass carried out at quite low tem¬ 
peratures (200-300 °C), atmospheric pressure and in an inert environ¬ 
ment - is currently considered as a key process in facilitating raw 
biomass (i.e., grass, woody, crop or processing industry residues) market 
growing for energy applications, due to its potential to convert any organic 
material into a high-energy-density solid fuel with superior properties in 
view of transport, handling, storage and in many major thermochemical 
end-use applications (e.g., co-firing in coal-fired power stations and 
gasification-based production of gaseous biofuels or bio-chemicals) [2], 

Olive pomace is the solid and wet by-product generated by olive-oil 
extraction industries. Among other possibilities of utilization (e.g. ani¬ 
mal feed, production of fertilizer, composting, etc.), the thermal treat¬ 
ment is appealing because it is an intensive and effective method for 
converting the olive pomace into energy [3], avoiding release of pollut¬ 
ants and odors as well as uncontrolled fermentation [4], Several studies 
[3,4] have been conducted on energy generation from the olive pomace; 
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however, no one includes torrefaction as pretreatment, which is able to 
produce an intermediate and stable material and does not suffer high 
moisture content, low energy density, hygroscopic behavior and limited 
storage life. 

With this background, the primary aim of the present study was to 
propose a set of kinetic parameters to be used as an input to future 
computer-based simulations of the olive pomace torrefaction, which 
can be used to guide laboratory experiments as well as to generate 
and evaluate alternative designs [5-7] for renewable energy projects. 

Non-isothermal thermogravimetric analysis (TGA) was selected as 
the experimental technique to obtain kinetic data due to its ability to 
overcome one of the major problems associated with standard isother¬ 
mal experiments, i.e. the warm-up time. Actually, during the initial non- 
isothermal time of an isothermal TGA test the sample typically un¬ 
dergoes some transformations that are likely to affect the results of 
the subsequent kinetic analysis [8], In addition, when a too slow heating 
rate is applied, the weight loss taking places during the warm-up phase 
is not negligible, which complicates the interpretation and deduction of 
kinetic data [9], 

Even though traditional model fitting approaches, including single 
or multi-step reaction mechanisms and/or single or multiple pseudo¬ 
component models [9-13], have found to date the almost exclusive ap¬ 
plication in biomass thermochemical conversion kinetic analysis [14], a 
“model-free” approach based on integral isoconversional methods was 
adopted in this paper. This choice was motivated by the fact that the 
model fitting procedures, when applied to non-isothermal data, pro¬ 
duce Arrhenius parameters strongly depending on the adopted reaction 
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model. As a consequence, they are so uncertain and variable that they 
cannot be meaningfully used for predicting the behavior of a substance 
over the range of experimental temperatures [8], Conversely, model- 
free techniques do not require any assumption of reaction models, 
thus avoiding possible errors associated with such a choice. In addition, 
yielding kinetic parameters as a function of either conversion 
(isoconversional analysis) or temperature (non-parametric kinetics) 
[15], these methods make it possible to study the overall kinetics 
throughout the whole range of experimental conversions and tempera¬ 
tures also for complex solid state processes, like biomass thermal de¬ 
composition, which usually involve several steps with different 
activation energies. 

Among the large number of differential and integral isoconversional 
methods to date developed, two integral multiple temperature- 
programmed methods, i.e., the linear Ozawa-Flynn-Wall (OFW) and 
the nonlinear Vyazovkin incremental methods, were selected to be 
comparatively applied in this work for the kinetic analysis of olive pom¬ 
ace torrefaction from non-isothermal TGA data. It is well-known [8, 
16-18] that, in a comparison between the two methods, the OFW pro¬ 
cedure is computationally simpler, but less accurate. 

In addition to proposing a set of kinetic parameters for a future con¬ 
ceptual process design of olive pomace torrefaction, this work aims to 
provide a further insight into the application of the model-free 
isoconversional methods, most commonly used for the kinetic analysis 
of thermally simulated process involving polymers [19] or other com¬ 
plex organic compounds [20], also in the field of biomass thermochem¬ 
ical conversion processes. 

2. Materials and methods 

2.1. Feedstock sampling and characterization 

The olive pomace used in this research work was suitably sampled 
from the accumulated residues of a three-phase olive mill located 
in Frasso Telesino (41° 9' 22 " N, 14° 31' 48" E) in Campania region, in 
November 2013. In order to preserve their original qualities and to pre¬ 
vent the microbial degradation before further analyses, all samples, just 
collected, were sealed in a plastic bag and stored at —20 °C Prior to use, 
representative fractions of this agro-industrial waste were milled using 
a knife mill to a particle size lower than 1 mm, afterwards, if required, 
oven-dried at 105 °C for about 12 h and then stored in a desiccator. 

Proximate (TG 701 LECO Thermogravimetric analyzer), ultimate 
(CHN 2000 and SC 144 DR LECO analyzer) and calorific (Parr 6200 
Calorimeter) analyses were performed in order to obtain the basic prop¬ 
erties as a fuel of this agro-industrial residue. The results are presented 
in Tablet. 

The indicative chemical composition of the olive pomace from a 
three-phase centrifugation process is also shown in Table 1, based on 
data reported in the literature [21], 

2.2. Experimental methods 

Thermogravimetric analyses were carried out in a TA Instruments 
analyzer Q600 SDT. Nitrogen was used as the purge gas at a flow rate 
of 100 ml/min to ensure an inert atmosphere and to prevent secondary 
reactions by volatiles produced during the solid thermal decomposition. 
During each run, approximately 10-15 mg of oven-dried samples, with 
a size of 0.5-1 mm, were placed in an aluminum oxide crucible inside 
the furnace chamber of the microbalance. Dynamic runs at five different 
heating rate programs, i.e., 2,5,10,20,40 °C/min, were carried out over 
the temperature range from room temperature to about 1000 °C to be 
used for the kinetic computations. 

In addition to this series of constant heating rate runs, two isother¬ 
mal experiments at 250 °C and 300 °C were also performed to check 
the reliability of the employed kinetic approach. During the isothermal 
run, the dried sample was first heated up to a temperature of 105 °C, 


Table 1 

Proximate, ultimate, calorific and chemical analyses of olive pomace. 


Moisture (wt.%, as received) 

Proximate analysis (wtX , dry basis) 
Volatile matter 
Fixed carbon 


Ultimate analysis (wt% dry basis) 


Ash 

O (by diff.) 

HHV (MJ/kg, dry basis) 

LHV (MJ/kg, dry basis) 

Chemical composition (wt.% dry basis) [21] 

Hemicelluloses 

Celluloses 

Lignin 

Fats and oils 


Olive pomace 
62.71 


52.58 

0.78 

7.23 

<0.01 

37.23 

21.88 

20.49 


15.91 

34.90 

20.51 

7.82 

6.89 


which was held for 10 min. Then the temperature was increased to 
the desired test value at a heating rate of 5 °C/min and kept constant 
for 7 h. Note that the aforesaid isothermal time excludes the warm-up 
phase, which globally lasted for about 1 h. Both thermogravimetric 
and differential temperature measurements were recorded simulta¬ 
neously during each test 

2.3. Non-isothermal kinetic analysis 

All the integral curves derived from dynamic heating rate experi¬ 
ments (TGA) were used in the kinetic analysis, after being converted 
to the plot of conversion (a) versus temperatures by Eq. (1) 


_ m(Tj)—m(T) 
m(Tj)—m(T f ) 


where m(T), m(Tj) and m(Tf) designate the instantaneous sample 
weight, the initial sample weight and the final sample weight, respec¬ 
tively. To provide the kinetic analysis with a standard reference, the ini¬ 
tial and the final sample weights were identified, for each run, at 150 °C 
and 450 °C, which correspond to the temperature range where olive 
pomace decomposition mainly occurs (see detailed discussion in 
Section 3.1 ). The motivation for this choice lies in the fact that, although 
the major objective of torrefaction is substantially the complete degra¬ 
dation of its hemicellulose content, it is unavoidable that such decom¬ 
position takes place without a simultaneous and partial degradation of 
the other biomass components, in particular the cellulose fraction [22], 

As mentioned in Introduction, the model-free kinetic approach 
based on two integral isoconversional methods, i.e., the Ozawa- 
Flynn-Wall (OFW) and the Vyazovkin ones, was employed in this 
study in order to evaluate the effective activation energy (E„) as a func¬ 
tion of the conversion degree (a). 

Although it is traditionally expected that a kinetic analysis produces 
a description of a process in terms of the so-called kinetic triplet [8], 
i.e., the reaction model and Arrhenius parameters (f(a), E and A), no 
further computations aimed at evaluating the pre-exponential factor 
and the reaction model were performed in this study. The reason is 
that the sole goal of the kinetic analysis carried out in this study was 
to predict the characteristic time needed to achieve a prefixed conver¬ 
sion degree under conditions typical of torrefaction that could be a 
valuable result for designing a torrefaction reactor. As shown by 
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Vyazovkin et al. [17], the knowledge of the E„ vs a dependence is by 
itself sufficient to do this. 

Following the pragmatic guidance from the ICTAC Kinetics Commit¬ 
tee [17], the reliability of the global activation energy obtained from the 
dynamic TGA runs, was checked by simulating isothermal profiles, ob¬ 
tained from tests not used in the kinetic analysis. To provide a standard 
reference for all the experimental variables dependent on the conver¬ 
sion degree by both the isothermal and non-isothermal TG, during 
this validation step, the aforesaid independent variable, which, by defi¬ 
nition, is a normalized form of the experimental weight loss, was re¬ 
defined as follows: 


It is worth noting that, although the use of approximation of the tem¬ 
perature integral instead of an accurate numerical method induces an 
error in the value of the activation energy, such approximations are typ¬ 
ically preferred since they lead to the simplest final plots (i.e., linear) 
yielding the activation energy very quickly [16], 

The application of the very crude Doyle's approximation [25] leads 
to the following linear Eq. (6) that is at the base of the popular linear 
Ozawa-Flynn-Wall method used in this work: 

,n ft- to (|^)- 5 - 3305 - , 052 ift; < 6 > 


_ m(ti)-m(t) 

~ m(t,) 

where m(t) and m(q) designate the instantaneous sample mass and the 
initial sample mass (i.e., sample mass at 150 °C), respectively. 

Simulation of the experimental curves was obtained by means of the 
following expression [12]: 



where t„* is the time to reach the conversion degree a * at an arbitrary 
temperature T 0 (i.e., 250 °C and 300 °C) under isothermal conditions, 
P = dT / dt is the constant heating rate, T is the absolute temperature, 
T a * is the temperature at which the conversion a* is experimentally 
achieved at the heating rate p, R is the gas constant and E„* is the acti¬ 
vation energy related to the extent of conversion, a*. 

The integral in Eq. (3) was solved based on the values of E a * obtained 
by the application of the OFW method and experimental T„* by dynam¬ 
ic thermal curves recorded at the same heating rate of the simulated iso¬ 
thermal ones (i.e., 5 °C/min). Specifically, in the present study, Eq. (3) 
was used to compute the time at which different conversion degrees 
were reached at 250 °C and 300 °C, under isothermal conditions. 

2.4. Ozawa-Flynn-Wall and Vyazovkin isoconversional methods for the de¬ 
termination of the apparent Ea dependence on a 

Model-free isoconversional methods rest upon the isoconversional 
principle, which states that, at a constant extent of conversion, the reac¬ 
tion rate is a function only of the temperature [19], 

Different integral isoconversional methods originate from the vari¬ 
ous ways of integrating the following single step kinetics Eq. (4): 


where T„j is the temperature at which the conversion a is experimen¬ 
tally achieved at the heating rate p, and the subscript i is a number 
representing experiments performed at different heating rates, Pi. 
Through the Eq. (6), this method allows to obtain the value of the appar¬ 
ent activation energy E a at a fixed conversion value a, from the slope of 
the straight line generated in the plot of log (Pi) versus 1/T^i, for a num¬ 
ber n of TGA experiments carried out at different heating rates (Pi). 

In this study, by determining the value of E H for each conversion de¬ 
gree in the range of a = 0.05-0.95 with a step of 0.05, the dependence 
of E„ on a was obtained. 

According to the more accurate nonlinear Vyazovkin method, for 
each given value of a, the apparent activation energy value, E„, can be 
determined by minimizing the following function ©(Ed), 


^) = EtrE^ 


EaJa,i)'Pj 


(7) 


where the temperature integrals are solved numerically, with a very 
time-consuming procedure. 

In this work, the standard Solver function, contained in the Microsoft 
Excel spreadsheet package, was used to obtain the value of E a minimiz¬ 
ing the function <t>(E H ). This solver applies the Generalized Reduced 
Gradient (GRG) nonlinear solving algorithm that is one of the most ro¬ 
bust nonlinear programming methods [26], For each run, in particular, 
after substituting experimental values of T a and p into the nonlinear 
Eq. (7), the value of each temperature integral was found by applying 
the trapezoidal integration rule with a temperature step of 10“ 2 K 
and assuming the constancy of E a only for a small integral conversion 
Aa ( = 0.05) according to the formula [19]: 

r(E a ,T a ) = J^exp^-^dT (8) 


W =AeXp (“^)' f( “ ) (4) 

where t is the time, A is the pre-exponential factor, E is the activation 
energy, T is the absolute temperature, R is the gas constant, a is the con¬ 
version degree of the process that can be determined from TGA runs as a 
fractional mass loss and f(a) is the reaction model, that is a function de¬ 
pending solely on a. 

For a constant heating rate program, where the temperature is 
raised at a constant rate p = dT / dt, the integration of Eq. (4), in partic¬ 
ular, yields: 


Note that the original method [27] uses the regular integration from 
0 to T h , as a result of which each value of E tt becomes averaged over the 
region 0 — a, and the whole E a dependence undergoes an undesirable 
flattening. However, the methods were subsequently modified to ade¬ 
quately account for a variation of E„ with a [28], To do this the regular 
integration from 0 to T„ was replaced with the integration over small 
temperature segments as in Eq. (8). 

3. Results and discussions 

3.1. Thermogravimetric analysis 


where g(a) is the integrated form of the reaction model f(ct) and I(E, T) is 
the temperature integral. Since the temperature integral in Eq. (5) has no 
analytical solution, it is typically solved by using approximations or nu¬ 
merical integration, each of which is at the base of a different integral 
isoconversional method. 


Mass loss (TG) and derivative mass loss (DTG) curves for the ther¬ 
mal decomposition of olive pomace at different heating rates (i.e., 2, 5, 
10,20,40 °C/min) are shown in Fig. 1. 

As expected, due to the presence of residual oil that inevitably and in a 
variable percentage, remains trapped within the matrix of this residue 
after the extraction process, the olive pomace does not show the typical 
thermogravimetric behavior of a standard lignocellulosic biomass (mainly 
composed of hemicellulose, cellulose and lignin), which is characterized 
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Temperature (°C) 

Fig. 1. Olive pomace (a) TG and (b) DTG curves at different heating rates (2, 5,10,20, 
40 °C/min) recorded in nitrogen atmosphere with a purge rate of 100 ml/min from ambi¬ 
ent temperature to 1000 °C 


by three decomposition phases, in addition to the first one due to the re¬ 
moval of water and water soluble components in the sample. In fact, as 
shown by the thermogram recorded with a heating rate of 2 °C/min 
(Fig. 1 ) and zoomed in greater detail for better readability in Fig. 2, five dis¬ 
tinct weight loss phases could be identified on the TG curve, related to the 
four downward peaks observed on the DTG curve, being the last three 
partially overlapping, and to the final gently sloping baseline. The first 
phase, corresponding to the reduction in mass (ft; 4.5 wt.% of the total 



mass loss) at temperature lower than 150 °C, can be attributed to the 
removal of water and light volatiles. The second phase, between 150 
and 275 °C, which accounts for approximately 23.4 wt.% of the total 
mass loss, can be attributed mainly to the thermal degradation of hemicel- 
lulose fraction, which is the most reactive part of a biomass and typically 
decomposes in the range of 220-280 °C [22], and partly to the beginning 
of cellulose decomposition. The third phase (275-395 °C) which accounts 
for approximately 40 wt.% of the total mass loss, in turn can be attributed 
mainly to the thermal degradation of cellulose, which typically degrades 
in a higher temperature range of 240-350 °C [22] and partly to the end 
of hemicellulose decomposition. In the fourth phase, over the temperature 
range from 395 to 450 °C, there was a negligible weight loss (ft; 6.9 wt.% of 
the total mass loss) that can be attributed to the thermal decomposition of 
the olive pomace oily fraction. Finally, the fifth phase located between 450 
and 1000 °C corresponds to the slow (i.e„ without a characteristic peak) 
weight loss of lignin whose degradation typically occurs over a wider tem¬ 
perature range [23] as a consequence of the different thermal stability of 
the various oxygen functional groups from its structure [24], This stage 
contributes to a weight loss of 25.2 wt.% of the total mass loss. 

The characteristics related to all decomposition phases such as the 
starting temperature, Tj, the ending temperature, T f , the weight loss, 
the maximum weight loss rate, W max , and the temperature where this 
rate occurred, T max , were exactly determined for all the studied heating 
rates (i.e., 2,10, 20,40 °C/min) and their values are given in Table 2. 

It is well-known that the heating rate affects the TG and DTG curves 
considerably [29], More specifically, it can change the curve shape and 
characteristic temperatures indicated in a curve. Fig. 1 shows the heating 
rate effects on the TG (Fig. 1 (a)) and DTG (Fig. 1 (b)) curves of olive pom¬ 
ace. It appears that higher heating rates gave globally a lower weight loss 
(or conversion). In addition, when the heating rate increases, TGA curves 
are shifted toward the right and the peaks of DTG curves are slightly 
shifted toward higher temperatures. Accordingly, all the characteristic 
temperatures in Table 2 were shifted to higher values with increasing 
heating rate. The findings observed in this study mirror those of the pre¬ 
vious investigations that have examined the same effects [30], 

Fig. 3 shows the weight loss (TG) curves and derivative of the mass 
loss curves (DTG) obtained from two isothermal runs performed at 
250 °C (a) and 300 °C (b), under inert atmosphere. It also includes the 
temperature versus time plots describing the temperature program 
used. At 250 °C over a test time of 475 min (Fig. 3(a)), the total weight 
loss of olive pomace sludge is 36.42 wt.%. The most evident and consid¬ 
erable weight loss (~49 wt.% of the total mass loss), however, occurs in a 
narrower time range (approximately 40-75 min) as highlighted by the 
spike in DTG, most likely as a result of the predominant decomposition 
of hemicellulose. 

Note that, despite the warm-up phase that is relatively short in this 
test due to the relatively low temperature (250 °C) adopted, the weight 
loss during this stage is not negligible, as it accounts for approximately 
30 wt.% of the total mass loss. This supports what has been said previ¬ 
ously about the limitations of a kinetic analysis based on isothermal 
data, with regard to the error introduced by the sample decomposition 
during the warm-up phase. As indeed observed in the present work, the 
weight loss during the non-isothermal period of an isothermal experi¬ 
ment may be considerable, despite a relatively short warm-up phase 
occurs, since it also depends on the specific reactivity of the sample. 

Moving from 250 °C to 300 °C and keeping the two runs of equal 
length, the total conversion of olive pomace increases by about 20.87% 
(Fig. 3(b)), about 50% of the total mass converted being lost during 
the warm-up phase. The onset of a second peak in the DTG curves 
displayed in Fig. 3(b), confirms that at more severe torrefaction condi¬ 
tions (300 °C), cellulose is more deeply degraded. 


Temperature (°C) 


32. Global kinetic analysis 


Fig. 2. Olive pomace TGA (••■) and DTG (—) curves at p = 2 °C/min recorded in nitrogen 
atmosphere with a purge rate of 100 ml/min. Vertical dashed-dotted lines delimit the 
five distinct weight loss phases. 


The OFW plots of ln(pj) versus 1 /T (l i for different values of conver¬ 
sion (a) are shown in Fig. 4. For each value of a, the apparent activation 
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energy was calculated from the slope of the plotted regression line. The 
results obtained and the correlation coefficient (R 2 ) are listed in Table 3 
together with the values of the apparent activation energy E„ obtained 
by using the nonlinear Vyazovkin procedure. 

By comparing the aforesaid values of E a (Table 3), it results that the 
error in the activation energy, as induced by the approximation of the 
temperature integral used to derive the simplest linear equation 
(i.e., Eq. (6)) on which the OFW method is based, may be considered 
reasonably small, being lower than the conventionally accepted 10% 
level in the activation energy [8], Hence, with regard to the olive pom¬ 
ace decomposition, this means that the linear procedure may be conve¬ 
niently used for the computation of reliable kinetic parameters. The 
greatly reduced computational load as compared to the more accurate 
nonlinear Vyazovkin incremental procedure just offsets a slightly higher 
error in the activation energy determination. 

Fig. 5 shows the dependence of the activation energy on the extent 
of olive mill sludge conversion, as computed by both the linear and non¬ 
linear methods, in a temperature range of interest for a torrefaction pro¬ 
cess. A small decrease in the E„ appears at low conversion degree 
(<0.15), followed by an almost constant (no multiple peaks and/or 
shoulders) value of E„ over the conversion range 0.2 < a < 0.70 and 
by a subsequent rapid rise of E H for a ranging between 0.75 and 0.95. 

It is worth noting that such trend reflects, in a way, the three main 
steps of thermal decomposition of olive mill sludge (Fig. 1(b)), falling 
in the range of 150-450 °C where the kinetic evaluation particularly fo¬ 
cused on. A comparison of the two trends, based on the values of the 
conversion degree, reveals that the first trend change observed in 
Fig. 5 at a = 0.15 corresponds roughly to the maximum of the second 
peak of the DTG curves in Fig. 1(b) as well as the second trend change, 
at a = 0.7, corresponds roughly to the maximum of the third peak of 
the DTG curves in Fig. 1(b). Consistently, the lower the value of activa¬ 
tion energy in Fig. 5 the higher the weight loss rate in Fig. 1(b). 

Anyway, the apparent activation energy as given by both the linear 
(i.e., ranging between 167.94 and 994.09 kjmol -1 ) and nonlinear 
(i.e., ranging between 162.59 and 1069.74 kj-mol -1 ) procedures global¬ 
ly shows the same increasing dependence upon the conversion degree 
in the range of 0.05-0.95. This trend can be qualitatively explained by 
the fact that the residual biomass becomes increasingly more reluctant 
to further decomposition in a N 2 atmosphere as the conversion 
progresses. 

Note, however, that the lower value of the correlation coefficient at 
a = 0.90-0.95 (Table 3) implies that the corresponding values of E a 
might be less reliable. 

As it is well-known, during the torrefaction process, a biomass 
typically loses 20-35 wt.% of its initial mass, mainly as a conse¬ 
quence of the complete degradation of hemicellulose and of the 
partial decomposition of cellulose. Accordingly, this means that 
over the range of interest for torrefaction (i.e., a < 0.55 based on 
Eq. (2)) in Fig. 5, the apparent activation energy calculated for the 
thermal decomposing of olive mill sludge may be regarded as 
roughly constant, amounting to about 200 kj mol -1 by averaging 
E a over a. This finding suggests that the torrefaction of olive pomace 
can be adequately described by a single-step model, which provides 
an adequate kinetic representation of a multi-step process having a 
single rate-limiting. Note that the averaging of E a is a procedure that 
may be justified only when a change in E re as a function of a is just a 
few percent units of the mean value, but not when such a change is 
comparable to the mean value [8], 

On the basis of the apparent activation energy dependence on the 
extent of decomposition, as obtained by the OFW method, curves 
representing the conversion of olive pomace at two different isothermal 
conditions (i.e., 250 °C and 300 °C) were obtained by the use of Eq. (3). 
To give an indication of the reliability of the kinetic approach employed, 
the results of the above prediction were then compared with data ob¬ 
tained by the isothermal TG analysis (Fig. 3), not included in the kinetic 
computations. 
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Fig. 3. TGA (-) and DTG (-) i 
b) 300 "C 



| 250 °C and 


It is worth observing that, according to the ICTAC Kinetics Commit¬ 
tee [17], the validation of the computed kinetic parameters by demon¬ 
strating that they can be used to satisfactorily predict (reproduce) 
experimental curves not included in computation rather than experi¬ 
mental curves used to derive the same parameters is a more rigorous 
approach. 



ent values of conversion. 


Table 3 

Activation energy for olive pomace decomposition by the Vyazovkin and Ozawa-Flynn- 
Wall methods. R 2 corresponds to the linear fitting in Fig. 4. 


a Vyazovkin Ozawa-Flynn-Wall method 
(—) method 


E„ E„ Plot equation R 2 Relative 

(kj mol -1 ) (kj mol -1 ) error 3 


05 

10 

15 

20 

25 

30 

35 

40 

45 

50 

55 


70 

75 


90 

95 


lnp = -29,463 T -1 
lnp = - 23,760 T -1 
lnp = -21,404 T -1 
lnp = -21,250 T -1 
lnp = -22,506 T -1 
lnp = -24,303 T -1 
lnp = - 26,050 T -1 
lnp = -27,005 T -1 
lnp = - 27,345 T -1 
lnp = -27,643 T -1 
lnp = -27,745 T -1 
lnp = -27,812 T -1 
lnp = -27,631 T -1 
lnp = -27,762 T -1 
lnp = -27,762 T -1 
lnp = -35,140 T -1 
lnp = -49,093 T -1 
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Fig. 5. Dependence of the activation energy on the conversion degree as determined by the 
isoconversional Vyazovkin and Ozawa-Flynn-Wall methods. 

As shown in Fig. 6(a)-(b), predictions are in excellent agreement 
with results from isothermal measurements. It is recalled that the dot¬ 
ted lines in Fig. 6(a)-(b) were derived from the isothermal data 
shown in Fig. 3(a)-(b), respectively, whereas simulation results 


(i.e., circle symbols) were calculated on the basis of the dynamic TG 
curves (Fig. 1(a)), via the OFW procedure. 

Note that, since the above simulations were performed by using 
Eq. (3) that allows to compute the time at which a given conversion is 
reached at an arbitrary temperature under isothermal conditions, the 
comparison was focused only on the isothermal region of the thermo¬ 
grams at T = 250 °C (Fig. 3(a)) and at T = 300 °C (Fig. 3(b)) and, 
more specifically, limited to the narrow time range (t < 150 min) 
where the most rapid weight loss occurs. 

4. Conclusion 

Decomposition kinetics of olive pomace in the temperature range of 
interest for torrefaction were investigated by the non-isothermal TG at 
different heating rates (i.e., 2, 5,10, 20,40 K/min). The dependence of 
the apparent activation energy on the conversion degree was deter¬ 
mined by means of a model-free approach based on two integral 
isoconversional methods. The reliability of results was checked by sim¬ 
ulating two isothermal experimental profiles (i.e., at 250 and 300 °C) 
that were not used in the kinetic analysis. The predictions were in excel¬ 
lent agreement with the experimental data. This means that the 
torrefaction of olive pomace can be adequately described by a single- 
step model. The findings of the present study might open perspectives 
for the implementation of a computational model aimed at determining 
optimal torrefaction conditions with respect to the olive pomace. 
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Fig. 6. Comparison of predicted data and experimental curve for olive pomace decompo¬ 
sition at: a) 250 °C and b) 300 "C. 
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